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Evolution of the cluster temperature function is extremely sensitive to the 
mean matter density of the universe. Current measurements based on cluster 
temperature surveys indicate that Qm ~ 0.3 with a la statistical error ~ 0.1, 
but the systematic errors in this method are of comparable size. Many more 
high- z cluster temperatures will be arriving from Chandra and XMM in the near 
future. In preparation for future cluster temperature surveys, this paper analyses 
the cluster mass-temperature relation, with the intention of identifying and 
^ ' reducing the systematic errors it introduces into measurements of cosmological 

parameters. We show that the usual derivation of this relation from spherical 
top-hat collapse is physically inconsistent and propose a more realistic derivation 
based on a hierarchical merging model that more faithfully reflects the gradual 
ceasing of cluster evolution in a low-f^M universe. We also analyze the effects of 
current systematic uncertainties in the M v - n —Tx relation and show that they 
introduce a systematic uncertainty of ~ 0.1 in the best-fitting Q M . Future 
improvements in the accuracy of the M vir —Tx relation will most likely come 
from comparisons of predicted cluster temperature functions with temperature 
functions derived directly from large-scale structure simulations. 
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1. Introduction 

Surveys of distant clusters of galaxies are now realizing their promise as cosmological 
indicators (e.g., Henry 1997; Eke et al. 1998; Borgani et al. 1999a; Donahue & Voit 1999). 
Because clusters are the largest virialized objects in the universe, and the latest objects to 
form in hierarchical models of structure formation, their rate of evolution is quite sensitive 
to cosmological parameters. However, because cluster masses are difficult to measure 
directly, a surrogate for cluster mass is usually used when comparing cluster observations 
to structure-formation models. In the X-ray regime, the simplest cluster observables are 
X-ray luminosity {Lx) and emissivity- weighted temperature (Tx). Temperature is more 
directly related to a cluster's mass, but luminosity can also be mapped to mass via an 
Lx — Tx relation. In either case, the cruical link between models and observations is the 
mass-temperature relation. 

Recent analyses of high-redshift cluster temperature functions from the Einstein 
Extended Medium-Sensitivity Survey (EMSS) have shown that the matter density of the 
universe probably lies in the range 0.2 < Qm < 0.7 (Henry 1997; Donahue et al. 1998; 
Bahcall & Fan 1998; Eke et al. 1998; Donahue & Voit 1999; but see Blanchard & Bartlett 
1998; Viana & Liddle 1999). Even though these conclusions are based on rather few clusters, 
the systematic errors in these measurements of f2 M are comparable to the statistical errors 
(Donahue & Voit 1999). With the flood of cluster temperature measurements expected over 
the next few years from Chandra and XMM, we will have the opportunity to measure Qm 
much more precisely. If we are to take full advantage of these measurements, we will need 
to reduce the systematic errors that currently exist in the modeling. Here we concentrate 
on the uncertainties in the mass-temperature relation itself. 

This paper analyzes the role of the mass-temperature relation in characterizing cluster 
temperature evolution. Section 2 outlines the formalism used to decribe the evolution of the 
cluster mass function. Section 3 investigates the physics underlying the mass-temperature 
relation, showing that the standard derivation of this relation from spherical top-hat 
collapse is flawed and suggesting a new context for understanding the physical effects that 
govern this relation. Section 4 discusses how to normalize the mass-temperature relation 
and evaluates how severely uncertainties in this normalization affect measurements of 
cosmological parameters. Section 5 summarizes the paper. 
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2. Mass-Function Evolution 

Numerical simulations have shown that the Press-Schechter formalism (Press & 
Schechter 1974) and its various extensions (e.g., Lacey & Cole 1993) characterize 
gravitationally driven structure formation with surprising fidelity, particularly on cluster 
scales (e.g., Lacey & Cole 1994; Eke, Cole, & Frenk 1996; Bryan & Norman 1998; Borgani 
et al. 1999a). Because of its successes, this formalism is often used to relate observed cluster 
temperature and luminosity functions to cosmological models (e.g., Henry & Arnaud 1991; 
Eke et al. 1996; Borgani et al. 1999a). While some recent large-scale simulations have 
shown that the Press-Schechter formula might be somewhat less successful at predicting the 
number density of the most massive clusters (e.g., Governato et al. 1999), we will assume in 
this paper that it is an exact description of the cluster mass function, because here we are 
more concerned with analyzing systematic problems with the mass-temperature relation. 

This section briefly outlines the formalism we will employ for expressing the cluster 
mass function. We derive expressions for mass-function evolution in both open and flat 
universes, and we assess the effects of a cosmological constant on cluster evolution. Many 
of these results have been derived elsewhere; we compile them here as background for 
subsequent sections. 

2.1. Press-Schechter Formalism 

The Press-Schechter formalism for structure formation and its extensions describe how 
virialized objects grow from a field of initial density perturbations. One defines <5(x, t; M) to 
be the local fractional overdensity of the universe smoothed on mass scale M and centered 
on comoving point x at time t. While these fluctuation amplitudes are linear, they grow in 
proportion to the function D(t), which depends on Qm and Q^, and their rms amplitude on 
scale M can be expressed as a(M)D(t)/D(t ). Ultimately, some of these fluctuations grow 
non-linear, and they are assumed to virialize when their amplitudes, extrapolated from 
the linear regime according to D(t), exceed some critical threshold for virialization 5 c (t), 
which also depends on Qm and Q^. One can then trace the merger history of a mass parcel 
beginning at location x from time t\ to the present by keeping track of the largest M for 
which 5(x,t 1 ;M)D(t)/D(t 1 ) > 5 c (t). 

Assuming that the perturbations are Gaussian, we can assess the number 
density of virialized objects with mass > M by evaluating the quantity 
u c (M,t) = [5 c (t)/a(M)][D(to)/D(t)], which is the critical virialization threshold in 
units of the characteristic fluctuation amplitude. The probability that a given mass parcel 
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is part of a virialized structure of mass > M is then equal to erfc[i/ c (M, t)/y/2\, where 
erfc(g) is the complementary error function. Thus, the overall mass density in virialized 
objects exceeding mass M is 




p(> M) = poerfcl^l = ^ / e" x * dx 



oo 2 



7T Jv c /V2 

where po is the mean mass density of the universe. Differentiating this expression with 
respect to M and dividing the result by M yields the familiar Press-Schechter formula for 
the comoving differential number density dn of virialized objects within mass interval dM: 



dM y ' ' \tx) M 2 



d In a 



d\a.M 



v. 



: (M,t)exp[-^ 2 (M,t)/2] , (2) 



where p cr) o = 3Hq/8itG represents the present-day critical mass density. 

The evolution of dn/dM depends solely on u c (M,t). At present, we have 
u c (M,to) = 5 c (to)/a(M), and in principle, we can determine a(M) by fitting equation (fj) 
to the current distribution of cluster masses. The function a(M) can be approximated by a 
power law with index a = (n + 3)/6 on cluster scales, so that a(M) = a 8 (M/M s )~ a with 
M 8 = (H$tt M /2G)(8 h- 1 Mpc) 3 = 6.0 x 10 14 ^m^ 1 M . Holding a(M) fixed, we can project 
the current cluster mass distribution backward in time as long as we know the functions 
8 c (t) and D(t). More specific expressions for a(M) describe how n changes with the mass 
scale in CDM-like models, but the analytical simplicity of the power-law form better serves 
the illustrative purposes of this paper. 



2.2. Cluster Evolution and 

Massive cluster evolution is very sensitive to Qm because the number density of large 
clusters depends exponentially on v 2 which is ^> 1. Slight differences in the rate at which 
v c evolves therefore translate into large differences in cluster evolution (e.g., Oukbir & 
Blanchard 1992; Eke et al. 1996; Viana & Liddle 1996). To illustrate how dramatic these 
differences can be, we will briefly outline the case of cluster evolution in an open universe 
with no cosmo logical constant. 

Following Lacey & Cole's (1993) treatment of perturbation growth when < 1 and 
Q\ = (see also the Appendix), we have u c (M,t) = uj(t)/a(M) with 



w(t) = 5 c (t)D(t )/D(t) = |r>(fo)[l + (V*) 2/3 ] , 



(3) 
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where tn = 7iQ M /H (l — fi M ) 3 / 2 . Normalizing the present value of v c at some fiducial mass 
scale M , we can write 



Because the cluster number density at mass scale M obeys dn/dM oc z/ c exp(— u%/2), we 
can use this equation to gauge how rapidly clusters at this mass scale evolve. 

Figure [l] illustrates how sensitively the number-density evolution of massive clusters 
(M « 5 x 1O 14 /i" 1 M ) depends on ft M . For the purposes of this illustration, we have 
adopted the <r 8 fitting formulae of Eke et al. (1996) and have assumed that n = —1.5 
on this mass scale. In this example, the number density of clusters in a universe with 
f^M — 0.8 grows by about five orders of magnitude from z = 1 to the present, a rapid rate of 
evolution that contrasts sharply with the single order of magnitude expected in a universe 
with Qm = 0.2. Changing the perturbation spectrum within observationally allowed bounds 
changes the quantitative predictions somewhat, but the qualitative conclusion remains the 
same: the number density of massive clusters evolves much more rapidly for Qm ~ 1 than 
for fl M 1. 



Clusters evolve slightly more rapidly in a flat, tt\ > universe than they do in 
an open, Q A = universe with the same value of f2 M (e.g., Eke et al. 1996). When a 
cosmological constant is operating, the universe's density remains close to critical later 
in time, promoting perturbation growth at lower redshifts. However, cluster evolution is 
considerably less sensitive to Q\ than it is to Qm- 

In order to characterize cluster evolution in a flat universe with Q\ > we require 
expressions for D(t) and S c (t), which are derived in the Appendix. From these expressions 
we can construct the threshold function 



where £ c (£), defined in the the Appendix, is proportional to the specific energy of a 
perturbation that collapses at time t. Normalizing v c as before at the mass scale M , we 
can write 





2.3. 



Cluster Evolution and A 




(5) 




(6) 



Plugging this expression into equation (fj) then yields the desired formulae for cluster 
evolution. 
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The dotted lines in Figure [I] show how the number density of clusters at mass scale M 
evolves when > 0, for the <7g normalization of Eke et al. (1996) and n = —1.5. Note 
that the rate of cluster evolution is quite insensitive to In this particular case, the 
best-fitting a% for a flat universe is slightly higher than that for an open universe, which 
almost compensates for the slightly more rapid rate of evolution owing to A. In general, the 
best-fitting Qm in an open universe is ~ 0.1 higher than in a flat universe (e.g., Donahue & 
Voit 1999). 

3. Theoretical Mass-Temperature Relations 

The Press-Schechter formalism conveniently describes the rate at which virialized 
objects of mass M accumulate in the universe. If we could observe cluster masses directly, 
then comparisions between Press-Schechter predictions and observed cluster evolution 
would be simple. Several types of observables, such as X-ray temperature, cluster velocity 
dispersion, and weak lensing, are related to cluster masses, but linking these quantities with 
the proper Press-Schechter M values requires careful attention. 

This section focuses on the relation between cluster mass and cluster temperature, 
the crucial relationship for linking X-ray observations of clusters to models of structure 
formation. We first outline the observational evidence for a well-behaved mass-temperature 
relationship. Then we analyze the standard derivation of the mass-temperature relation, 
which is based on collapse of a spherical top-hat perturbation to an isothermal sphere. This 
derivation yields a relation similar to the observed relation, but it fails to conserve energy, 
indicating that it omits important physical effects. In an effort to understand this relation 
more deeply, we present a model for cluster formation, drawn from the merging-halo 
formalism of Lacey & Cole (1993), which accounts for the fact that massive clusters accrete 
matter quasi-contiuously. Analyzing clusters in this context enables us to identify the 
physical effects that make up for the lack of energy conservation. The primary advantage of 
the continuous formation model is that it more naturally reproduces the late-time evolution 
of clusters, and the section concludes by comparing predictions of cluster temperature 
evolution drawn from the continuous formation model with those from the spherical top-hat 
model. 
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3.1. 



Observational Evidence 



Simple scaling arguments suggest that the X-ray temperatures of clusters (Tx) should 
be directly related to their masses. One way to define a cluster's mass is to specify a 
characteristic radius r& within which the mean density is A times the critical density p cr , 
so that Ma = 47rr A p cr A/3. If all cluster potentials share the same density distribution, 
p(r/r A ), and the X-ray gas is isothermal, then T x oc M A /r A oc M A 3 . Numerical simulations 
of cluster formation demonstrate that this scaling ought to be remarkably tight, with a 
scatter of only 15 - 20% (Evrard, Metzler, & Navarro 1996; Bryan & Norman 1998). These 
simulations also provide normalizations for the mass-temperature relation that can be 
compared with actual clusters. 

A recent observational investigation of the cluster mass-temperature relation at z ^ 0.1 
by Horner, Mushotsky, & Scharf (1999) supports the results of the simulations. They 
show that cluster masses derived from velocity dispersions (Girardi et al. 1998) agree well 
with those inferred from ASCA temperatures (Fukazawa 1997), using the scaling law from 
Evrard et al. (1996) at A = 200: 



The scatter in the observed mass-temperature relation is ~ 30%, but it decreases by a factor 
of 2 for the clusters with the highest numbers of measured galaxy redshifts, suggesting that 
the scatter intrinsic to the mass-temperature relation is probably quite small. However, 
some concerns remain: a handful of outliers deviate from the standard relation by up to 
50% and the mass normalization one finds from hydrostatic modeling of a subset of these 
clusters is 40% smaller (Horner et al. 1999). 

Similar comparisions at higher redshifts are more difficult, but the available data 
indicate that the mass-temperature relation remains well-behaved. Hjorth, Oukbir, & van 
Kampen (1998) have compared masses derived from gravitational lensing analyses for 8 
clusters at 0.17 < z < 0.54 with the X-ray temperatures of these clusters. Their best-fit 
mass-temperature relation agrees with the Evrard et al. (1996) scaling law within the 
observational errors, and they conclude that this scaling law can be used to measure masses 
to within 27%. For clusters at even higher redshifts (0.53 < z < 0.83), Donahue et al. 
(1999) show that the observed relation between X-ray temperatures and cluster velocity 
dispersions remains consistent with the \ow-z relation. 




3/2 



(7) 
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3.2. Virial Mass and the Late-Formation Approximation 

The seemingly good behavior of the cluster mass-temperature relation is fortunate for 
those who wish to study cluster evolution with X-ray telescopes, but care must be taken 
when relating these temperature-derived masses to the virial masses demanded by the 
Press-Schechter formalism. From the simulations and observations, we know that the mass 
within a specified density contrast is straightforwardly related to temperature. However, 
the density contrast A vir corresponding to the virial radius depends in general on Qm 
and fl\. Thus, in order to characterize cluster evolution properly, we need to know how 
A vir (t; fl M , ^a) changes with time. 

The usual approach to defining a cluster's virial mass is to approximate cluster 
formation with the evolution of a spherical top-hat perturbation (e.g., Peebles 1993). Such 
a perturbation formally collapses to the origin at a particular moment (t c ) which is taken to 
be the moment of virialization. The virialization time thus equals twice the time required 
for the perturbation to reach its turnaround radius (r ta ). A naive application of the virial 
theorem, assuming that the perturbation is cold at maximum expansion, dictates that 
the cluster's final potential energy ought to be twice its potential energy at turnaround. 
Hence, the cluster's virial radius is assumed to be half its turnaround radius (r vir = r ta /2). 
According to this prescription, A vir is a well-defined function of cosmic time and the 
parameters f^M and (Lacey & Cole 1993; Kitayama & Suto 1996; Oukbir & Blanchard 
1997). In the case of £l\ = 0, this function can be concisely expressed as A vir = 8-rc 2 /(Ht) 2 , 
where H is the Hubble constant at time t. If we additionally assume that each cluster we 
see at a given redshift z has just reached the moment of virialization, an assumption known 
as the late- formation approximation, then M vir oc T^/ 2 Pc^ 2 ^iJ 2 ■ 

In a critical Qu = 1 universe, the late-formation approximation is valid because 
massive clusters develop rapidly at all redshifts; the effective moment of virialization is 
always close to the moment of observation. However, in a universe with Qm < 1, cluster 
formation is currently shutting down, and one must account for differences between the 
moment of virialization and the moment of observation. This problem grows most severe 
at late times in a Qm <C 1 universe, because the quantity p cr A vir as determined via the 
late-formation approximation declines indefinitely. The M vir associated with a given T x 
therefore rises steadily, even though cluster evolution has essentially stopped. This spurious 
late-time evolution of the M vir —T x relation is an undesirable artifact of the late-formation 
approximation. 

One approach to solving this problem is to account explicitly for the difference between 
the moment of virialization and the moment of observation in the context of a merging-halo 
formalism for cluster growth (Viana & Liddle 1996; Kitayama & Suto 1996). Another 
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equivalent but mathematically simpler approach, which § |3.3| describes in detail, is to 
consider how the M vir —T x relation should evolve in a population of clusters that gradually 
accrete their matter over an extended period of time, a more realistic scenario for the 
growth of very massive clusters (Voit & Donahue 1998). 

Calculating the normalization of the M vir — Tx relation under the late-formation 
approximation is also somewhat problematic. Because a virialized cluster's potential is 
approximately isothermal, one would like to approximate it with a singular isothermal 
sphere, truncated at radius r vir , within which the mean density is p cr A vir . The one- 
dimensional velocity dipersion within such a potential is cr 2 D = GM/2r\ lv (Binney & 
Tremaine 1987), which leads to the following relation between virial mass and temperature: 



kTv — ■ 



2(3 



47T 



1.38 keV)/T V^Mi/X/ 3 



n 



1/3 



(8) 



where (3 = nm v (j\ Q jkTx and jj,m p = 1 x 10~ 24 g is the mean mass per gas particle. 

Comparisions of the M Y - U —Tx relation in equation (H) with the masses and temperatures 
of simulated clusters indicate that (3~ x ~ 0.8 — 1 (e.g., Bryan & Norman 1998). This 
nearness of (3 to unity appears to validate the assumptions governing the derivation of 
equation (||), but the approximate agreement between this equation and the simulations 
turns out to be something of a coincidence. The total energy of a collapsing spherical 
top-hat perturbation is — 3GM 2 /5r ta . After the collapsed perturbation virializes into an 
isothermal sphere, a naive application of the virial theorem that disregards boundary 
effects would place the total kinetic energy of the system at 3GM 2 /5r ta , corresponding to 
af D = 2GM/5r ta . The virial radius of the relaxed system would then be 5r ta /4, a factor 
of 2.5 larger than assumed in the derivation of equation (|8|), and its temperature would 
correspondingly be 2.5 times lower. 

In fact, truncation of a virialized system at some r vir implies the existence of a 
confining pressure, unaccounted for in the top-hat collapse model, that alters the usual 
virial relationship between potential and kinetic energy (e.g., Carlberg, Yee, & Ellingson 
1997). In the case of a singular isothermal sphere, the total kinetic energy is three times 
the absolute value of the total energy. Thus, energy- conserving collapse of a spherical 
top-hat perturbation into a pressure-truncated singular isothermal sphere should yield 
<r 2 D = 6GM/5r ta , implying r vir = 5r ta /12 (e.g., Shapiro, Iliev, & Raga 1999). This result is 
close to the naive assumption that r vir = r ta /2, but it is valid only if a confining pressure is 
applied at the virial radius. 
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3.3. Continuously Forming Clusters 

The inconsistencies in the top-hat, late- formation derivation of the M wir —Tx relation 
outlined above indicate, perhaps unsurprisingly, that the top-hat collapse model excludes 
important physical effects that contribute to the normalization of the relation. In particular, 
the top-hat model accounts for neither the energy and mass accumulation during the early 
stages of cluster formation nor the confining effects of matter that continues to fall in, both 
of which significantly increase the temperature associated with a given mass. This section 
shows how these missing effects can be addressed in the context of a simple model in which 
massive clusters are allowed to form gradually, rather than instantaneously. 

In hierarchical models for structure formation, the growth of the largest clusters is 
quasi-continuous. The most massive clusters are so rare that they almost never merge with 
another cluster of similar size (e.g., Lacey & Cole 1993). Rather, they grow by continually 
accumulating much smaller virialized objects. In the notation of § their masses grow like 
M oc uj- 3 /( n+ V (Lacey & Cole 1993; Voit & Donahue 1998). Because each bit of infalling 
matter carries with it a specific energy e, we can compute the virial energy —E of the 
cluster by integrating E = — J edM. The cluster temperature itself is proportional to E/M, 
so this integral also leads to a relation between virial mass and temperature. 

Voit & Donahue (1998) treat the case of continuous cluster growth when < 1 and 
= 0, finding that M oc x~ 3m ^ 5 , where x = 1 + (tn/t) 2 ^ 3 and m = 5/(n + 3). Here we 
extend that calculation to include the constant of proportionality between energy and mass. 
Drawing on the Appendix, we express the specific energy of infalling matter at time t as 

1 /27rGM\ 2/3 



Thus, we obtain 

E 3 m /2nG\ 2/3 



j M 2 ' 3 ( 



'tn\ 2/3 + l_ 
t J m 



(10) 



M 10 m- 1 V tn 

In the limit of large m, which corresponds to the late-formation approximation, this 
expression reduces to 

- = ~-e(t) , (11) 
M 5 v ; 

which is identical to the E/M ratio for a spherical top-hat perturbation of mass M that 

virializes at time t. A similar procedure yields the mass-temperature relation in a flat 

> universe. From the Appendix, we have e(t) oc M 2 / 3 £ c (t) and u(t) oc — £ c (i), giving 

E 3 m , , 
M 5 m — 1 
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which again reduces to — 3e(i)/5 in the limit of large m. 

Two factors in equation ( |TUD drive E/M higher than the late- formation value. The 
171/(171 — 1) factor, also present in equation (|l|), accounts for the effects of early infall; 
continuous cluster formation tends to create hotter clusters than top-hat formation because 
more of the mass is assembled early, at a higher mean density. For values of n typical 
of cluster scales (—2 ^ n ^ —1), this factor ranges from 1.2 to 1.7. The 1/m term in 



the bracketed factor of equation (|1(J) accounts for the cessation of cluster formation when 
t 3> tn. At late times in an open universe, E/M should remain constant, but in the 
late-formation approximation E/M falls indefinitely because the fiducial density scale never 
stops dropping. 

Relating E/M to temperature requires an expression for the relationship between the 
total virial energy —E and the total kinetic energy Ek- When an external pressure P 
confines the boundary of a spherically symmetric virialized system, the appropriate form of 
the virial theorem can be written 

E K = E + 4irPrl iv . (13) 
If we take the velocity dispersion to be isothermal (am = const.), then P = p(r vir )o" 2 D , and 

E K = _ P E , (14) 
p - 2p(r vir ) 

where p is the mean density within the virial radius. In this formulation, the ratio Ek/E 
depends on the shape of the potential within r vir . If the local density is negligible at r vir , 
then the confining pressure is effectively zero and Ek = E. If the potential strictly obeys 
p oc r~ 2 , then Ek = 3E. 

Because we wish to derive an approximate mass-temperature relation for comparison 
with observations and simulations giving the temperature and mass of a cluster within r '200, 
let us compute Ek/ E for the "universal" density profile of Navarro, Frenk, & White (1997), 
truncated at r 20 o, under the assumption that am is constant. The ratio of mean density to 
local density within r 2 oo is then 



P M + c) 



2 



3- ^ln 



p(r 200) c 2 



;i + c )- 



l + c 



(15) 



where c is a parameter that quantifies the concentration of matter toward the cluster's 
center. Simulations by Eke, Navarro, & Frenk (1998) show that c m 4 — 6.5 for clusters 
at < z < 1 in a — 0.3, — 0.7 universe, and simulations of cluster formation by 
Navarro et al. (1997) show that the most massive clusters in critical universes exhibit 
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similar levels of concentration. Values of c ~ 5 lead to density contrast factors ~ 4 and 
E K /E « 2. 

An intriguing alternative approach by Shapiro et al. (1999) investigates the post- 
collapse structures of clusters by seeking the minimum-energy solution among a family 
of non-singular truncated isothermal spheres. The minimum-energy solution turns out to 
closely resemble the self-similar spherical infall solution of Bertschinger (1985). Because 
the truncation radius of this minimum-energy isothermal sphere is nearly equal to the 
accretion-shock radius of the infall solution, Shapiro et al. (1999) suggest that continual 
infall naturally maintains the confining pressure on the virialized isothermal sphere. In this 
model, the density contrast factor at the truncation radius is 3.73, and E^/E = 2.17. 

Taken together, the effects of early accretion and pressure confinement make up for the 
lack of energy conservation in the top-hat, late-formation derivation of equation In the 
early-time limit (t <C t$-j), equation ( |i~0D yields the following mass-temperature relation: 



kT 



x 



m Ek 



5 m — 1 E 



GM 2 l*iim p /4vr 



2(3 V 3 



1/3 

p cr A vir ) . (16) 



When n ~ —2 and E^/E rs 2, the prefactor in brackets is close to unity, making this 
expression nearly identical to equation (H). 

The lesson here is that the assumptions underlying equation are physically unsound. 
The approximate agreement between the M V1T —Tx normalization derived via the top-hat 
collapse model and those derived from simulations and observations is largely coincidental. 



As long as f2 M is not very small, the time-dependent factors in equations @ and (|T0D do 
not differ by a large factor. However, because equation fli~0|) more faithfully reflects the 
behavior of cluster formation in all the appropriate limits, we prefer to base the M Y - 1T —Tx 
relation on the continuous-formation model. 



3.4. Late-Formation vs. Continuous-Formation 

The M vir —T x relations in equation (||), derived using the late-formation approximation, 
and equation (|10D, derived using the continuous- formation approximation, differ in both 
normalization and time-dependent behavior. The following section will discuss the 
importance of properly normalizing the M vir —Tx relation. Here we wish to examine how 
differences in time- dependence alone translate into different predictions for cluster evolution. 
In order to isolate the time-dependent behavior, we can identically normalize both M vir —Tx 
relations to equation (0) at z — and compare the resulting cluster temperature functions. 
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Figure |2| shows the evolution of the temperature function for 8 keV clusters, given the 
erg normalization of Eke et al. (1996). Because the M wir —Tx relation evolves less strongly in 
the continuous-formation case, the rise in cluster temperature at a given mass as z increases 
does not compensate as fully for the drop in the number of clusters at that mass. The 
evolution of the temperature function at a given value of Qu is therefore stronger in the 
continuous-formation case (see also Viana & Liddle 1999). Correspondingly, the best-fitting 
Qm to a given observed amount of cluster evolution will be lower. In this particular case, 
the difference amounts to ~ 0.1 in Qm for a best-fitting Qm ~ 0.3. 

Because the statistical errors in Qm derived from cluster temperature evolution 
are also ~ 0.1, this discrepancy between the late- formation and continuous-formation 
approximations will need to be resolved if we are to take full advantage of the cluster 
temperature measurements expected from Chandra and XMM. The best way to proceed 
will be to test how well these temperature-function predictions represent the results of 
large-scale structure simulations. However, cluster temperature functions will have to 
be extracted directly from the simulated data, without resorting to a mass-temperature 
conversion step, presumably by using the cluster-particle velocity dispersion as a surrogate 
for cluster temperature. 

4. Normalization of the M vir — Tx Relation 

Both the top-hat and continuous-formation derivations of the M y \ T —Tx relation given 
in the previous section have holes which must be plugged with knowledge gained from 
simulations. In fact, any such spherically symmetric representation glosses over aspects of 
cluster formation that are inherently three-dimensional. Thus, it seems wise to normalize 
these analytical expressions to the results of numerical simulations. This procedure appears 
simple enough, but one must bear in mind that the normalization depends on and that 
simulations have been done only for a few particular values of Qm- Here we explain how we 
choose to normalize the M wlx —Tx relation then investigate the consequences of an offset in 
the normalization. 



4.1. Normalizing to Simulations 

Because of the good agreement between the observational compilation of Horner et al. 
(1999) and the simulations of Evrard et al. (1996), we would like to normalize the M viT —Tx 
relation accordingly. Applying the time- dependence factors derived for continuously forming 
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clusters to the empirical mass-temperature relation in equation (0) thus gives 



2/3 

kT x = (8.0 keV) 



M~ V" (t n /t) 2/3 + l/m 



10 15 h' 1 Mr. 



(V*o) 2/3 + l/m 



(17) 



and 

kT x = (8.0 keV) 
for open and flat universes, respectively. 



M \ 2/3 Ut) 



10 15 h~ l M & J Ut ) 



Figure ^ compares this empirical normalization with the normalizations of the 
M vir — Tx relations derived by Eke et al. (1996) and Voit & Donahue (1998) at z — 0, 
assuming fl^ = 0. The dashed line indicates the normalization given in equation (0), 
which is presumed to be independent of Qm- The curve labeled "late formation" shows 
the normalization of equation (Kfi, derived from the top-hat, late- formation model. This 
normalization drops steadily with decreasing fl M because the density contrast factor A vir 
grows smaller as Qu declines. Clusters modelled in this way are therefore less compact 
and cooler than one would expect from the critical density alone. When Qm = 1, this 
normalization is only 4% below the empirical value, but if Qm — 0.2, it lies 20% below this 
value, corresponding to mass discrepancies of 6% and 30%, respectively. 

The behavior of the normalizations derived from continuous-formation models is more 
complicated. Voit & Donahue (1998) normalized these relations to the Eke et al. (1996) 
relation at Qm = 1 to simplify comparisons. For Om 1, they are lower than at f^M = 1 
for the same reason as in the late-formation model. However, if Qu *C 1, cluster formation 
happened long before z ~ 0, when the universe was considerably denser. Clusters are 
therefore denser and hotter than one would expect from the current critical density. As a 
result, the temperature normalization of the n = —2 case deviates by less than 10% over 
the range 0.2 < Qm < 1-0. In the n = —1 case the normalization is actually 18% higher 
than the empirical value (27% in mass) at Qm — 0.2. 

This comparison illustrates why the procedure of normalizing the M vir — Tx relation 
to simulations is imperfect. In general, we expect this normalization to vary with in a 
way that depends on n. Given that we have simulations for only a handful of cosmo logical 
models, how do we unambiguously normalize these relations? Furthermore, different choices 
for extrapolating this normalization to other values of can lead to normalizations that 
differ by as much as 40% in temperature (60% in mass) at a given value of £1m- Because 
of these uncertainties in the normalization of the M vir —Tx relation, it is important to 
understand how offsets in the normalization affect cosmological parameters derived from 
the cluster temperature function. 
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4.2. Normalization Offset and a 8 

The M vir —Tx relation is invoked twice in the usual derivation of as from the low-redshift 
temperature function. In both instances, an overestimate of the mass associated with a 
given temperature drives the best-fitting value of a§ higher. For example, a 50% offset in 
the mass normalization changes as by about 15%. Systematic uncertainties in the M^—Tx 
relation therefore lead to systematic uncertainties in cr 8 , limiting the usefulness of the 
temperature function as a tool to measure the perturbation amplitude. 

The first place the M v ; r — Tx relation enters is in the conversion of the theoretical 
cluster mass function dn/dM in equation (0) to the cluster temperature function 

z/ c [M(T,t),t]exp{-z/ c 2 [M(T,t),t)]/2} . (19) 

An M v i r — Tx relation that overestimates M(T) by a fractional amount 5m will underpredict 
the density dn/dT by the same fractional amount. This source of error drives the best-fitting 
value of v c lower by a fractional amount 5 V ~ —5m/ {v 2 ~ !)• 

Once v c {T) has been derived over a given range of temperatures, one can determine 
<t(T) = 5 c (t ) /p c {T). These a values will be too high by a fraction ps 5m /(y 2 ~ 1) if there 
is a normalization offset. Conversion of cr(T) to cr(M) contributes another term to the 
systematic error budget. If M(T) is overestimated, the mismapping of temperature to mass 
inflates as by a fractional amount ~ cxSm- 

As an example of these effects, consider the derivation of as from the abundance of 
> 5 keV clusters. Markevitch (1998) finds that the number density of such clusters at 
z ~ is ~ 7.0 x 10~ 7 h 3 Mpc -3 . According to the mass-temperature relation in equation 
(0), the mass of a 5 keV cluster is 5.0 x 10 14 hr 1 M & , leading to an overall mass density 
p(> 5keV) ~ 2.3 x 10 -32 h 2 gem -3 in such objects. Plugging this value into equation ([[]) 
and solving for v c yields z/ c (5keV) ~ 3.2 for Qm = 1 and z/ c (5keV) ~ 2.8 for Qm = 0.3, 
numbers that are consistent with more rigorous fits to cluster temperature data using a 
similar mass-temperature relation (Donahue & Voit 1999). 

Conversion of these v c values to as values depends on the shape of the initial perturbation 
spectrum and the value of the virialization threshold 5 C . For the purposes of this analysis, 
we will assume that 5 C is known perfectly, implying that <r(5.0 x 10 14 /i _1 M Q ) w 0.53 for 
f2 M = 1-0 and a(5.0 x 10 14 /i~ 1 M o ) w 0.59 for f2 M = 0.3. Extrapolating along the mass 
spectrum assuming, for example, that n = —1.5 then leads to as ~ 0.51 for = 1-0 and 
a 8 w 0.76 for Q M = 0.3. 

Now let us inflate the normalization of the mass-temperature relation by 50% in mass. 



dn (T ) = 3 /2\ l ' 2 fWcr.o 
dT [ ' ' 2 W TM(T,t) 
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The overall mass density in objects > 5 keV rises to p(> 5keV) xs 3.5 x 1CT 32 h 2 gem -3 , 
yielding z/ c (5keV) ps 3.1 for f2 M = 1 and z/ c (5keV) ps 2.7 for Q M = 0.3. Because 5 keV 
corresponds to 7.5 x 10 14 /i _1 M Q under the alternative normalization, we now have 
(j(7.5 x 10 14 /i" 1 M ) ps 0.54 for fi M = 1.0 and cr(7.5 x 10 14 h~ x M ) ps 0.61 for fi M = 0.3. 
Extrapolation to as, again assuming n = —1.5, gives a$ ~ 0.58 for Qm = 1-0 and a 8 ps 0.87 
for fl M = 0.3. Note that these systematic changes, of order 15%, exceed the typically 
quoted measurement errors for a 8 at a fixed value of Q M . 

Additional uncertainty in the derived value of as can arise from uncertainty in the 
slope n of the perturbation spectrum. For example, the n ps —2.3 slope derived from the 
Markevitch cluster sample (Markevitch 1998; Donahue & Voit 1999) leads to a considerably 
lower derived value of a 8 when Q M « 1. In the case of f2 M = 0.3, a measurement of 
cr(5.0 x lO 14 /i -1 M ) ~ 0.59 extrapolates to a 8 ~ 0.66, corresponding to clusters of 
temperature ~ 2.5 keV, below the temperature limit of the sample. Blanchard et al. 
(1999) have recently argued that the Markevitch sample is incomplete at low temperatures, 
implying n > —2.3. That is probably why the best-fitting a 8 values of Donahue & Voit 
(1999), based primarily on the Markevitch clusters, seem unusually low. (Errors on a 8 
quoted in that paper refer only to the statistical errors in a 8 at the best-fitting value of JIm-) 
Ideally, one would like to measure a 8 from the number density of clusters at temperatures 
corresponding to the appropriate mass scale, but to do this, one first needs a reliable 
M viT —T x relation, in addition to a well constrained value of f2 M . 

The upshot of this analysis is that a 8 values derived from the cluster temperature 
function contain systematic errors that depend on the mass-temperature relation. These 
systematic errors are currently comparable to the measurement errors. Until the M vir —T x 
relation is better understood, a 8 values derived from the cluster temperature function will 
have to be treated with caution. Conversely, predictions of cluster temperature functions 
that invoke a 8 values derived from other kinds of data will also contain systematic errors 
owing to normalization uncertainties in the M vir —Tx relation. 



4.3. Normalization Offset and Q M 

The systematic problems discussed above in relating a 8 to dn/dT need not lead to 
unwarranted pessimism about deriving Qm from the evolution of dn/dT. The key quantity 
in establishing the rate of cluster temperature evolution is not a 8 , but rather u c (T, t ), whose 
systematic errors are considerably smaller, ~ 5% instead of ~ 15%. In order to evaluate 
the impact of a M vir — Tx normalization offset on predictions of temperature evolution, 
we will first analyze the case of = 1, then consider how a normalization offset affects 
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measurements of 



M • 



When Q M = 1, the time- dependent parts of M(T,t) and v(T,t) simplify to 
M oc (1 + z)~ 3 / 2 and v c oc (1 + z)( 2 ~ 3a ^ 2 . The amount of cluster evolution at a fixed 
temperature T can therefore be written as 



*) = ^7^7 = (1 + *) (5_3a)/2 exp{-^^[(l + ^) (2 " 3a) - 1]} • (20) 



A fractional overestimate 5m of cluster masses thus leads to an underestimate of the amount 
of evolution by a factor ~ exp{5jvf[(l + z) 2 ~ 3a — 1]}. If the overestimate of cluster masses 
is 50%, this factor amounts to a 20% underestimate of evolution at z — 0.3 and a 70% 
underestimate at z — 0.8 for n = —1.5. 

These uncertainties are relatively minor compared to the expected amount of evolution. 
For example, if n = —1.5 and v c {T,z = 0) ~ 3.2, we expect C(T,z = 0.3) ~ 0.3 and 
C(T,z = 0.8) ~ 0.03. Larger values of u c , characteristic of hotter clusters, lead to even 
more evolution. Because these evolution predictions are over an order of magnitude larger 
than the systematic errors, M vir —Tx normalization discrepancies do not seriously affect the 
conclusion that cluster temperature evolution rules out Qm = 1, particularly when clusters 
at z > 0.5 are included. 

Partially because of the potentially significant uncertainty in cr 8 , certain authors have 
been cautious about the conclusions that can be drawn about from cluster temperature 
evolution (Colafranceso, Mazzotta, & Vittorio 1997; Viana & Liddle 1999; Borgani et 
al. 1999b). However, maximum likelihood methods of determining Qm that compare the 
cluster temperature function at z ~ directly with the cluster temperature function at 
higher redshifts (e.g., Henry 1997; Donahue & Voit 1999) can obtain stronger constraints 
on Qyi because they are differential measurements in which much of the uncertainty in 
as and 5 C cancels. Figure [| shows the cluster evolution predictions that result when a 
representative range of a% values is considered. Here we allow 0.5 < <J$0^ 7 ~ ' 1 M < 0.6. 
At the z 0.3 redshift of the Henry (1997) clusters the prediction of the high-og, = 1 
model is only a factor of two lower than the low-ag, Qm — 0.5 model, underscoring the need 
to identify systematic sources of uncertainty in <t 8 before deriving evolutionary predictions 
for clusters from it. However, Figure |5| paints a somewhat rosier picture. Here we allow 
v c o = ^ c (5keV, z = 0) to span a range that corresponds to a factor of two range in the 
mass normalization at 5 keV, or equivalently, a factor of two range in the number density 
of 5 keV clusters at z — 0. The resulting systematic uncertainty in the best-fitting Q M is 
£0.1. 
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4.4. Toward Greater Precision 

When using Press-Schechter methods to model cluster evolution, one should always keep 
in mind that they are useful because they efficiently approximate numerical simulations. 
Our confidence in these methods is rooted in the fact that they reproduce the mass function 
of simulated clusters reasonably well. Less work has been done on comparisons of simulated 
cluster temperature functions with temperature functions derived from Press-Schechter 
mass functions using an M v - 1T —Tx relation (e.g., Bryan & Norman 1998, Pen 1998). 

Given the ambiguities surrounding the M v - n —Tx relation and the very definition of a 
cluster's mass, the most robust way to model the evolution of dn/dT would seem to be 
with a version of the Press-Schechter formalism that describes the temperature function 
directly. In such a scheme, one would separate the function v c (T,t) that determines the 
evolution of dn/dT into a temperature-dependent part u c (T, 0) and a time-dependent part 
g(t;flM,^A)- In principle, g(t; Qm> ^a) could be derived from a grid of simulated cluster 
temperature functions. However, creating such a simulation set would be very expensive. 
The simulation volume would need to be extremely large to obtain adequate statistics on 
rare, high-temperature clusters. 

Maximum likelihood fits to cluster temperature surveys are essentially seeking the 
best-fitting u c (T,z) = h , c (T,0)g(z;QM,^A)- Figure |] shows that this technique is fairly 
robust with respect to systematic uncertainties in the M wir —Tx normalization. Statistical 
uncertainties in the normalization and slope of v c (T, 0) are handled naturally by the 
maximum likelihood method. Insofar as the dependence of z/ c (T, z) on cosmological 
parameters is accurate, this technique currently has the potential to deliver values of Qm 
that are accurate to ^ 0.1. However, it remains to be seen how accurately our assumed 
forms for g(z; Qm, ^a) reproduce the results of large-scale clustering simulations. 

5. Summary 

X-ray surveys of distant clusters are placing increasingly more stringent constraints on 
f2fvi- The lack of extreme evolution in the cluster temperature function strongly indicates 
that Qm < 1- One of the crucial ingredients in placing such constraints on Qm is the 
M v ; r — Tx relation that converts cluster temperatures to cluster masses, enabling us to 
relate X-ray temperature surveys to theoretical models for cluster formation. If we are to 
extract accurate values of f2 M from the larger cluster temperature surveys expected from 
Chandra and ASCA, we need to ensure that our M v ; r — Tx relation faithfully describes 
cluster evolution when coupled with Press-Schechter analysis. To that end, this paper has 
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analyzed our current understanding of the cluster mass-temperature relation in an effort to 
identify the systematic errors it introduces into measurements of cosmological parameters. 

We find that the usual derivation of the M vir —T x relation, which assumes that 
clusters form by spherical top-hat collapse and that we are observing them immediately 
after they formed, is physically inconsistent. The rough agreement between the M vir —T x 
normalization derived in this way and the normalization determined from numerical models 
of clusters is therefore somewhat coincidental. To obtain the proper normalization, one 
needs to account both for the fact that much of a cluster's mass accreted well before the 
moment we are observing it and for the non-zero density at r2oo, which requires a surface 
pressure term to be included in the virial theorem. 

Because of these shortcomings of the spherical top-hat picture, we advocate a more 
realistic scenario for deriving the M v - 1T —Tx relation in which clusters form quasi-contiuously. 
An expression for the M w - 1T —Tx relation can be derived in the context of hierarchical 
merging, but the normalization of this relation depends on the concentration parameter c of 
the cluster, which must be obtained from simulations. The primary advantage of this form 
for the M vir — Tx relation is that, unlike the spherical top-hat model, it properly reproduces 
the cessation of cluster evolution at late times if < 1- 

Given the systematic uncertainties in setting the normalization of the M vir —T x 
relation, we have investigated their impact on the derivation of cosmological parameters 
from the cluster temperature function. Because two applications of the M vir — Tx relation 
are needed to extract as, this parameter is particularly susceptible to uncertainties in the 
mass-temperature normalization: a mass-normalization uncertainty of 50% leads to a 15% 
uncertainty in <7g. However, only a single application of the M vir — Tx relation is needed to 
extract Qm, making it less vulnerable to normalization uncertainties. The systematic error 
in Qm owing to uncertainties in the M vir — Tx relation is ~ 0.1. 

Improvements in our understanding of the M vir —Tx relation await comparisions of 
theoretically-derived cluster temperature functions with structure-formation simulations 
large enough to contain many hot clusters. In essence, the dn/dT expression derived from 
the M v i r —Tx relation and the Press-Schechter mass function is no more than an elaborate 
fitting formula for representing the results of simulations. The ideal mass-temperature 
relation will therefore be the one that reproduces simulated temperature functions or 
velocity-dispersion functions most accurately. Even better would be a fitting formula that 
gives dn/dT directly without passing through murky intermediate steps involving ill-defined 
cluster masses. Until we have resolved these systematic uncertainties in deriving dn/dT, 
our constraints on Qm from cluster temperatures will not grow appreciably tighter. 
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APPENDIX 

CONTINUOUS CLUSTER GROWTH IN fl M < 1 UNIVERSES 

In the spirit of Gunn & Gott (1972), we can idealize continuous cluster growth as 
occurring through the sequential collapse and virialization of an infinite series of concentric 
shells, each obeying the equation of motion 

» GM AR 

where R is the radius of the shell encompassing mass M, and A is the cosmological constant. 
The specific energy of the matter in the shell is then 

R 2 GM AR 2 

which remains constant until the shell virializes. If the shell ever reaches the critical radius 
R C r = (3GM/A) 1 / 3 , cosmic repulsion dominates gravity from then on, and the shell never 
collapses. 

We can cast the equation of motion for the shell into dimensionless form by defining 
x = R/R C r, 9 = A l 'H, and f = ei£ 2 A _1 , so that 



dx 



2 „ x 2 ^' 2 



(3) 



A particular shell will collapse if x^ + 6£x + 2 = for some x in the range < x < 1. 
Solving this cubic equation, we find 

-o = 2 3 / 2 |e| 1 / 2 cos(|-|) , (4) 

where is defined by cosc^ = — (S^l 3 )" 1 / 2 with 7r/2 < < it. The shell therefore reaches 
its maximum radius R ci Xq(£) at a time A~ 1 / 2 ^ (^), where 



[MO x 1 ' 2 dx 
9 o(0 = / —z rm ■ ( 5 ) 
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Because of the symmetry of the motion, the shell collapses to the origin at 
t c (£) = 2A -1 / 2 # (O> which we will take to be the time of virialization. 

The overall scale factor of the universe obeys a similar equation of motion, for which 
the specific energy of the background matter is 

•2 tt2 tt2 

e b = y-^W-^ A a 2 (6) 
= M(i_q m _q a ) . ( 7 ) 

The universe changes from decelerating to accelerating when a = a cr = (fi M /2^ A ) 1/3 , and 
with the definition w = a/a cr , the dimensionless equation of motion for the background 
expansion becomes 



dw 

1e 



2 , w 2 
— + 2& + — 

6W 6 



1/2 

(8) 



where = e;,a cr 2 A 1 . 



Taking advantage of the formal similarities between these equations of motion, we can 
derive an expression for perturbation growth in an open universe. In the linear regime, 
Sp/p = —35w/w, where 5w = x — w, and at any moment in time we have 



Jo 



y l/2 dy r w y 1/2 dy 



At early times, when x 3 + 2> |6£x| and w 3 + 2> |6^k;|, we obtain 



(9) 



5 P aft ^ (^ 3 + 2) 1/2 r y 3/2 dy 
and at the earliest times (w C 1), this expression reduces to 

- = |(6-^>- (11) 
p 5 

Perturbations at early times grow like 1/(1 + z), as expected, and their amplitudes are 
proportional to the specific energy difference between the perturbation and the background. 
Note that if the universe is flat, the background specific energy vanishes, and the 
perturbation amplitude within a shell is directly proportional to the shell's specific energy. 

In the case of a vanishing cosmological constant, we can relate the perturbation 
amplitudes explicitly to their collapse times t c . When the cosmological constant is very 
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small, we have |£| 3> 1 and t c « (7r/3\/2)|£| 3 / 2 A x / 2 , so a shell that collapses and virializes 
at time t c carries with it a specific energy 



1 /2ttGM\ 2/3 



€= 2 



{—) • (12) 
If we define = 7rQ M / H (l — fi M — ^a) 3 ^ 2 , then at very early times, 



5p 3(12tt) 2 / 3 / M 2/3 



10 VW 



l + (W*c) 2/3 ' • (13) 



According to Lacey & Cole (1993), the growth rate for linear perturbations in this limit is 
D(t) ~ [(127r) 2/3 /10](t/tn) 2/ ' 3 ) and we retrieve their expression for S c (t): 

6 e (t) = lD(t)[l + (t a /t) 2 '*] . (14) 
Time therefore enters the Press-Schechter formula via the parameter v c oc 5 c (t)/D(t) oc 

[i + (W0 2/3 l- 

In the case of a flat universe with Q M < 1, we have 

j = -HD(t) (15) 

where 

nm m mi (u> 3 + 2)V 2 fW y W dy 

m = DHt)] = w 3/ 2 J, ¥ + 2) 3/2 > ( 16 ) 

in agreement with Eke et al. (1996). Inverting the function £ c (£) yields the function £ c (£) 
giving the specific energy e(t) = -R 2 r A£ c (t) of a shell that collapses to the origin at time t. 
The collapse threshold becomes 

S c (t) = -H c (t)D(t) , (17) 
and time enters the Press-Schechter formula via v c oc 8 c (t)/D(t) oc £ c (i). 
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Fig. 1. 



Expected evolution of massive clusters at mass scale M for different values of Qu- 
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Fig. 2. — Expected evolution in the comoving number density of 8 keV clusters for the late- 
formation approximation and the continuous formation approximation. If Q M ps 0.2, the 
best-fitting Qm for the late-formation and continuous- formation cases can differ by ~ 0.1. 
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Fig. 3. — Dependence of the M Y - U —Tx normalization on Qm- The dotted line represents the 
temperature of a cluster containing 10 15 h^ 1 M cluster within the radius r 2 oo, as predicted 
by the simulations of Evrard et al. (1996). The solid line labeled "late formation" shows the 
temperatures predicted by M vir —T x relation of Eke et al. (1996), and the other two solid 
lines represent the M vir —T x relations from Voit & Donahue (1998) for n = — 1 and —2. 
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Fig. 4. — Expected evolution in the comoving number density of 8 keV clusters for 
0.5 < agf^M 47-01 M < 0-6- Because of uncertainties in the M v - 1T —Tx relation, as is not 
tightly contrained by the cluster temperature function. If these uncertainties are taken at 
face value, they result in a considerable spread in predictions for cluster evolution. 
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Fig. 5. — Expected evolution in the comoving number density of 8 keV clusters over the range 
of v c values corresponding to the observed low-redshift cluster temperature function at 5 keV. 
Despite the uncertainties in the M vir —T x relation, the parameter z/ c0 = z/ c (5keV,,2 = 0) is 
fairly well constrained. A 50% change in the mass-normalization, or correspondingly, a 50% 
change in the comoving number density of 5 keV clusters, changes the best fitting v c by 
~ 5%. The resulting spread in predictions for cluster evolution is not as severe as one might 
have guessed from the uncertainties in <r 8 . 



